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Abstract 

We study discretization effects in a mixed-action lattice theory with domain-wall valence quarks 
and Asqtad-improved staggered sea quarks. At the level of the chiral effective Lagrangian, dis- 
cretization effects in the mixed-action theory give rise to two new parameters as compared to the 
lowest order Lagrangian for rooted staggered fermions - the residual quark mass, m res , and the 
mixed valence-sea meson mass-splitting, A m j x . We find that m res , which parameterizes explicit 
chiral symmetry breaking in the mixed-action theory, is approximately a quarter the size of our 
lightest valence quark mass on our coarser lattice spacing, and of comparable size to that of simu- 
lations by the RBC and UKQCD Collaborations. We also find that the size of A m i x is comparable 
to the size of the smallest of the staggered meson taste-splittings measured by the MILC Collab- 
oration. Because lattice artifacts are different in the valence and sea sectors of the mixed-action 
theory, they give rise to unitarity-violating effects that disappear in the continuum limit, some of 
which should be described by mixed-action chiral perturbation theory (MAXPT). Such effects are 
expected to be mild for many quantities of interest, but are expected to be significant in the case 
of the isovector scalar (ao) correlator. Specifically, once the parameters m res , A m ; x , and two others 
that can be determined from the light pseudoscalar meson spectrum are known, the two-particle 
intermediate state "bubble" contribution to the scalar correlator is completely predicted within 
MAXPT. We find that the behavior of the scalar meson correlator is quantitatively consistent with 
the MAXPT prediction; this supports the claim that MAXPT describes the dominant unitarity- 
violating effects in the mixed-action theory and can therefore be used to remove lattice artifacts 
and recover physical quantities. 
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I. INTRODUCTION 



Advances in both computers and algorithms now allow realistic nonperturbative calcu- 
lations of hadron masses and matrix elements from QCD first principles. In particular, the 
results of lattice simulations with three flavors of improved staggered dynamical quarks are 
in good agreement with experimental measurements for a wide range of light meson, heavy- 
light meson, and heavy-heavy meson quantities These successes, which include both 
"postdictions" such as the pion and kaon decay constants and predictions such as the 
shape of the D — > Ktv form factor [3j, indicate that many of the systematic uncertainties 
associated with lattice QCD calculations are now under control. Therefore lattice QCD can 
now be used to calculate reliably many weak matrix elements that cannot be determined 
experimentally but are important for flavor physics phenomenology. These matrix elements 
are needed to interpret the results of experiments such as those at the Tevatron, B-factories, 
and, in the near future, the LHC, and for precise tests of the Standard Model in the quark 
flavor sector. 

This paper is the first in a series leading to the calculation of the kaon bag parameter (4], 
Bk, with all sources of systematic error under control using the mixed-action lattice QCD 
simulation scheme pioneered by the LHP Collaboration 5] . Bk parameterizes the hadronic 
part of neutral kaon mixing, and, when combined with an experimental measurement of 
ex, constrains the apex of the Cabibbo-Kobayashi-Maskawa unitarity triangle. Because 
€k has been well-measured experimentally p], the dominant uncertainty in this constraint 
is that of Bk- It is likely that whatever new physics exists has additional CP- violating 
phases which will manifest themselves as inconsistencies between measurements that are 
predicted to be identical within the Standard Model. Thus a precise lattice determination 
of Bk helps constrain physics beyond the Standard Model and is an important goal of the 
phenomenology community. The mixed-action method, which we describe in detail in the 
following section, employs domain-wall valence quarks and improved staggered sea quarks, 
and is particularly well-suited for the numerical calculation of weak matrix elements such 
as B K . 

In order to understand and demonstrate control over the various sources of systematic 
error that enter the numerical calculation of Bk, we have divided this project into smaller 
pieces that are both interesting by themselves, and, when viewed as a whole, will lend 
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credibility to our Bk calculation. Because discretization errors will be one of the primary 
sources of uncertainty in Bk, in this first paper we study their effects in the mixed- action 
theory. In particular, we quantify the size of the explicit chiral symmetry breaking param- 
eter, m res , and the mixed-action parameter, A mix . We also test the range of applicability of 
mixed-action Chiral Perturbation Theory (MAXPT) by analyzing numerical mixed-action 
data for the isovector scalar correlator, in which unitarity-violating discretization effects are 
expected to be substantial for our choice of simulation parameters. In the second paper we 
will present results for the pseudoscalar decay constants, f n and fx- Because f n is well- 
known experimentally, it provides a test of the mixed-action lattice methodology, especially 
the combined chiral-continuum extrapolation using MAXPT. The ratio fx/ fir then allows a 



Only after we have 



model-independent determination of the CKM matrix element V u , 
tested the mixed-action method using the known quantity will we present results for Bk- 
These results will appear in a subsequent publication. 

In this work we study discretization effects in the mixed-action theory. We calculate the 
residual quark mass, m res . We obtain a value of m res ^2.7 MeV in the chiral limit on the 
a « 0.12 fm ("coarse") MILC lattices; m res is even smaller, approximately 1 MeV, on the 
a « 0.09 fm ("fine") MILC lattices. Given these values, we find that m res on the coarse 
lattices is a quarter the size of our smallest valence quark mass, and comparable to that of the 
RBC and UKQCD Collaborations on a similar lattice spacing [8] , indicating that the amount 
of chiral symmetry breaking is acceptably small in our mixed-action simulations. We also 
calculate the mixed valence-sea meson splitting, A m j x , on both the coarse and fine lattices. 
We find that a 2 A mix « (280MeV) 2 on the coarse lattices and a 2 A mix » (190 MeV) 2 on the 
fine lattices, which is about a factor of two smaller. Thus these results are consistent with 
the expected 0(a 2 ) scaling behavior of a 2 A m ; x , as well as the independent determination on 
the coarse lattices by Orginos and Walker-Loud 9J. We note that the numerical values of 
both m res and A mix are unique to this particular choice of mixed-action scheme, and would 
be different if, for example, one were to use highly-improved staggered quark (HISQ) gauge 
configurations for the sea sector [lol |. 

One might be concerned that the use of a mixed action could introduce new theoretical 
complications, and, in particular, sizeable unitarity-violating effects that are not described 
by mixed-action Chiral Perturbation Theory (MAXPT). Because the valence and sea quarks 
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in mixed-action simulations generically have different discretization effects, one cannot tune 
the valence-valence and sea-sea pseudoscalar mesons to remove all unitarity violations at 
nonzero lattice spacing. In other words, the mixed-action lattice theory is always partially 
quenched, and one cannot recover full QCD until the continuum limit. Thus, in order to 
extract physical quantities from mixed-action simulations, one must be able to account for 
all sources of lattice discretization effects (at a given order in MAXPT) and remove them. 
In quantities such as fx and Bk, unitarity- violating discretization effects are relatively 
mild, and only appear as meson mass-shifts inside chiral logarithms. One would like to 
demonstrate, however, that even when the unitarity-violation is more pronounced, it can 
still be described by MAXPT. 

In this work we show that MAXPT accurately describes the behavior of the isovector 
scalar 2-point correlation function. The ao correlator is particularly sensitive to unitarity- 
violating discretization effects in the mixed-action theory because it receives contributions 
from flavor-neutral two-meson intermediate states. The "bubble" contribution to the scalar 



correlator was calculated in MAXPT by Prelovsek llj. The expression contains only a few 
low-energy constants, all of which can be determined in fits to pseudoscalar meson mass 
data. Thus the size and shape of the bubble contribution to the ao correlator is completely 
predicted within MAXPT. We find that, for all valence-sea mass combinations on both the 
coarse and fine lattices, the predicted bubble contribution is quantitatively consistent with 
our mixed-action 2-point data. Thus we conclude that MAXPT describes the dominant 
unitarity-violating effects in the numerical simulations. In the case of the scalar, we cannot 
actually determine the ao mass with our existing data because, for most of our coarse lattice 
data, the size of the bubble contribution swamps the leading exponential contribution from 
the ao meson over almost the entire time range. In the case of Bk, however, unitarity- 
violating non-analytic discretization effects are predicted to be less than a percent of the 



continuum value of Bk over the relevant extrapolation range on the coarse MILC lattices [12 ]. 
Thus we can use MAXPT to remove these effects and recover precisely the continuum value 
of B K . 

This paper is organized as follows. In Section [Til we describe the actions and parameters 
used in our lattice calculations. We calculate the residual quark mass, m Tes , in Section lTOl and 
the mixed-meson splitting, A mix , in Section HVl Both of these quantities appear in MAXPT 
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expressions, and we use our value of A m ; x to analyze the behavior of the do correlator 
within the framework of MAXPT in Section |V] In Section [VI] we summarize our results and 
conclude. 



II. MIXED-ACTION LATTICE SIMULATIONS 



In this section we motivate the use of mixed-action lattice simulations for the determina- 
tion of hadronic observables, and, in particular, weak matrix elements. We then describe the 
details of the lattice actions and parameters used in our numerical mixed-action calculations. 



A. Theoretical advantages 

Any numerical lattice simulation in which the actions for the valence and sea quarks 
are different is a mixed-action calculation. In this work, however, we use the terminology 
"mixed- act ion" specifically to refer to the choice of actions used by the LHP Collaboration, 
i.e., a domain- wall valence quark action on top of a staggered sea. 

Staggered fermions are computationally cheap as compared to other standard discretiza- 
tions. Thus they allow the lightest dynamical quark masses and finest lattice spacings 



currently available in lattice simulations 13(. This numerical affordability comes at the 



cost, however, of additional theoretical complications. Staggered quarks come in four de- 
generate species called tastes [141]; consequently each flavor of staggered meson exists in 
sixteen tastes. Although the sixteen light pseudoscalar meson tastes are degenerate in the 
continuum, the symmetry that relates them is broken at nonzero lattice spacing, leading 



to discretization errors that are of 0(a 2 ) 15|. Because these errors are numerically sig- 



nificant at current lattice spacings, staggered lattice data must be fit to functional forms 
that include taste-breaking discretization effects calculated in staggered Chiral Perturbation 
Theory (SXPT) Q, Q, Q Q • The unphysical effects can then be subtracted to recover the 
desired physical quantity. 1 The staggered chiral effective Lagrangian contains a number of 



1 Although there is no rigorous proof, there is considerable supporting evidence that the rooting procedure 
used in numerical simulations with staggered quarks is correct. We therefore work under the plausible 
assumption that the continuum limit of the rooted staggered theory is QCD. Both numerical and theo- 



retical evidence supporting the rooting procedure are reviewed by Sharpe in Ref. 19] and by Kronfeld in 



Ref. 



20] 



6 



low-energy constants in addition to those of continuum XPT. In many quantities of interest 
such as masses and decay constants, these new parameters primarily enter next-to-leading 
order (NLO) chiral expressions in a straightforward way by producing an additive shift to the 



meson masses that appear inside chiral logarithms [161. 1171] . In other quantities (such as Bk) 
that are not protected by a partially-conserved current, however, taste-symmetry breaking 
leads to a more complicated operator mixing pattern that must either be accounted for with 
fully nonperturbative operator renormalization, which is difficult, or with the inclusion of 
many new operators in the SXPT calculation, which introduces many new undetermined 
coefficients in the joint chiral-continuum extrapolation 2l| . 

The calculation of weak matrix elements such as Bk with domain-wall quarks 



22 



23(, on 



the other hand, is theoretically simpler than that with staggered quarks because domain-wall 
quarks do not occur in multiple species. Furthermore, they retain exact chiral symmetry 
at nonzero lattice spacing up to exponentially small corrections, the size of which can be 



controlled by the length of the 5th dimension 



. Consequently, while the AS = 2 lattice op- 



erator still mixes with wrong- chirality operators, there are significantly fewer operators than 
in the staggered case, and nonperturbative renormalization can be applied more easily 24]. 
Lattice simulations with domain-wall fermions, however, are computationally more expen- 
sive than those with staggered fermions with comparable masses and lattice spacings 25]. 

The LHP mixed-action scheme retains the primary advantage of staggered lattice simu- 
lations while reducing the complications associated with taste-symmetry breaking. Thus it 
is particularly well-suited for the numerical calculation of Bk- The mixed- action method 
allows for light dynamical quark masses and fine lattice spacings, but at the computational 
cost of a quenched domain-wall lattice simulation. The chiral symmetry of the domain- 
wall valence quarks eliminates mixing with operators of other tastes and minimizes mixing 
with those of other chiralities, thereby making the nonperturbative renormalization more 
tractable. Furthermore, it makes NLO mixed-action XPT expressions, including that for 



Bk, largely continuum-like 12, 



26 



27 



28. 



29(]. 



B. Numerical details 



We now describe the details of our numerical mixed-action simulations. A summary of 
the ensembles and valence quark masses used in this work can be found in Table [B 
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TABLE I: Parameters of the MILC improved staggered gauge configurations and domain-wall 
valence quark propagators used in this work. Columns three and four show the nominal up/down 
(mi) and strange quark (m s ) masses in the sea, along with the corresponding pseudoscalar taste 
pion mass. Columns five and six list our partially-quenched valence quark masses, along with our 
lightest available domain- wall pion mass. 







sea sector 


valence sector 




a(fm) 


L 


ami/am s am^ 


am v 




0.09 


28 


0.0062/0.031 0.14789(18) 


0.0062, 0.0124, 0.0186, 0.046 


0.1201(14) 


0.09 


28 


0.0124/0.031 0.20635(18) 


0.0062, 0.0124, 0.0186, 0.046 


0.1213(15) 


0.12 


24 


0.005/0.05 0.15971(20) 


0.007, 0.02, 0.03, 0.05, 0.065 


0.1733(12) 


0.12 


20 


0.007/0.05 0.18891(20) 


0.01, 0.02, 0.03, 0.04, 0.05, 0.065 


0.1971(11) 


0.12 


20 


0.01/0.05 0.22447(17) 


0.01, 0.02, 0.03, 0.05, 0.065 


0.19931(94) 


0.12 


20 


0.02/0.05 0.31125(16) 


0.01, 0.03, 0.05, 0.065 


0.1938(21) 



We use the MILC lattices with 2+1 dynamical flavors of Asqtad-improved staggered 



fermions 



30, 



L4J ; a detailed description of the simulation parameters can be found in Refs 
3lJ. We have data on both the "coarse" (a w 0.12 fm) and "fine" (a » 0.09 fm) MILC 
ensembles, which have physical volumes ranging from approximately (2.5 - 3 fm) 3 . For each 
ensemble, the mass of the dynamical strange quark is close to its physical value, and the 
masses of the up and down quarks are degenerate. The ratios of the nominal up/down quark 
mass to the nominal strange quark mass in the sea sector are given in the the third column 
of Table E 

We generate domain-wall valence quark propagators 



22 



231 ] using the Chroma software 
, we HYP-smear the lattices 



system for lattice QCD [321 ] . Like the LHP Collaboration 
using the standard parameters given in Ref. 33| in order to reduce the size of explicit chiral 
symmetry breaking and proximity to the Aoki phase 34|. We also use the same values for 
the domain- wall height, M 5 =1.7, and extent of the fifth dimension, Lg=16. Using these 
parameters, we find that the residual quark mass is acceptably small; details of our m res 
calculation are given in Section IIHI 



We choose bare domain-wall quark masses in order to best suit the numerical calcula- 



S 



tion of Bk, although we are also using the resulting propagators for quantities such as the 
pseudoscalar decay constants. We therefore need a wide range of valence quark masses that 
allows us both to extrapolate our lattice data to the physical up/down quark mass and to 
interpolate to the physical strange quark mass. In order to keep finite-volume effects under 
control, however, we restrict the quantity m n L > 4. Based upon these constraints, we have 
chosen the domain- wall valence quark masses listed in the fifth column of Table [B Note 
that, unlike in other mixed staggered sea, domain- wall valence lattice simulations we have 
not made any attempt to tune the domain-wall pion mass in the valence sector to match the 
lightest staggered pion mass in the sea sector. This is because, although many possible tun- 
ings exist, none completely removes unitarity- violating effects from the mixed-action theory 
at nonzero lattice spacing; full QCD is only recovered in the continuum limit. We have, 
instead, generated many partially-quenched points, and will use the appropriate expressions 



in mixed-action Chiral Perturbation Theory 12|, |26|] to extrapolate our pseudoscalar decay 
constant and Bk data in future work. 

We have generated Coulomb gauge-fixed wall-source propagators with periodic and an- 
tiperiodic boundary conditions starting at timeslice zero for use in the calculation of the 
K° — K° matrix element. In order to make best use of the available computing time, we 
have only generated quark propagators on every fourth MILC gauge configuration; thus all 
of our data points are relatively uncorrelated. These wall-source propagators are used for 
our determination of the residual quark mass, m res , in Section IIHI In order to compute 
the parameter A mix , we formed mixed valence-sea mesons by tying the wall-source prop- 
agators together with Coulomb gauge-fixed (but not HYP-smeared) point source Asqtad 
valence quark propagators made with the MILC code 35]. The results for A mix are given 
in Section IIV1 We have also generated periodic-boundary-condition Landau gauge-fixed 
point-source propagators for nonperturbative operator renormalization using the method of 



Rome- Southampton 



36]. These propagators were used for our preliminary analysis of the 



isovector scalar 2-point correlation function in Ref. 37| because a point source is required for 
comparison to the MAXPT expression in Ref. It turned out, however, that relatively 
large statistical errors in the point-point scalar data made it difficult to resolve the small 
bubble contribution on the fine lattices. We therefore generated additional Coulomb gauge- 
fixed Z2 random-wall source propagators on a subset of our ensembles in order to reduce the 
size of the statistical errors in the scalar correlator data. The random-wall source gives the 
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same normalization as a point source (up to a factor of the spatial volume), but averages 
over all of the spatial points on the source timeslice. The analysis of the a® correlator with 
the improved data is presented in Section [V] 

In order to convert dimensionful quantities determined in our mixed-action lattice sim- 
ulations into physical units, we use the MILC Collaboration's recent determination of the 



scale r\ = 0.318(7) fm [13j. The Sommer scale-setting method [38J] has the advantage that 



the ratio r\/a can be determined precisely from a fit to the static quark potential [3l|, |39 |. 
By converting all of our data from lattice spacing units into r% units before performing 
any chiral fits, we account for slight differences in the value of the lattice spacing between 
ensembles. MILC obtains the absolute scale ri by first using the HPQCD Collaboration's 
calculation of a" 1 on the coarse and fine lattices from Upsilon spectroscopy 4^] to get r x on 
the coarse and fine lattices, and then extrapolating these values to the continuum. 



III. THE RESIDUAL QUARK MASS 

In this section we determine the value of the chiral symmetry breaking parameter, m res , 
in our mixed-action simulations. 

Domain-wall fermions respect flavor symmetry and have an approximate chiral symmetry 



22 



23 



411 ] . This is possible because of the introduction of an extra dimension; the Dirac 
operator for domain-wall fermions is a five- dimensional operator, with free boundary con- 
ditions in the fifth dimension. Light chiral fermion modes are exponentially bound to the 
four- dimensional surfaces at the ends of the fifth dimension. In the limit that the separation 
between the two domain walls, commonly denoted L s , is taken to infinity, the chiral sym- 
metry becomes exact. In practice the value of L s must be finite in numerical simulations, 
vpically between 10-20 lattice spacing units, and chiral symmetry is not exactly maintained 



As in the case of Wilson fermions, the leading chiral symmetry breaking terms are pro- 
portional to the lattice spacing a. However, for Wilson fermions this term is 0(1), whereas 
for domain- wall fermions this term is (9(10~ 3 ) for typical choices of L s . The only effect of 
this leading-order term is a shift of the effective bare quark mass by a small residual mass, 
m res . The generic, weak-coupling behavior of chiral symmetry breaking leads to the expec- 
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tation that m res depends on L s as aexp(— aL s ), for some constant a [42]. 2 Additional chiral 
symmetry breaking effects are suppressed by higher powers of the lattice spacing, on top of 
the exponential in L s suppression factor. Therefore the leading corrections to domain-wall 



fermions are 0(aexp(— aL s )) and 0(a 2 ) 44|. 



In practice, the residual mass is measured in lattice simulations using the ratio of the 



midpoint current to the axial current 45|, |46j : 



As t — > oo, this quantity should approach a constant which is then defined as the residual 
mass, that is, R(t) — > m res . The axial current in Eq. ([1]) is 

J 5 a (x) = L s - l)P L t a V(x, 0) + V(x, 0)P R t a V(x, L.-l), (2) 

and the midpoint current is 

Jsqfr) = L s/2 ~ l)P L t a ^(x, L./2) + L s /2)P R t a V(x, L s /2 - 1), (3) 

where $ is the domain-wall field, the flavor matrices t a are normalized by Tr(t a t b ) = 5 ab , 
andP flii = (l± 75 )/2. 

When defined as in Eq. (Q]), m res is constant up to discretization errors which give m res a 
dependence on the light quark masses (both valence and sea) and the lattice spacing. This 
light quark mass dependence has been observed in a number of earlier simulations. In our 
simulations, we follow the method of the LHP Collaboration by applying HYP smearing to 
the gauge fields before inverting the domain- wall propagators |5j. This smearing reduces 
the coupling of the domain-wall quarks to high momentum gluons which exchange chirality, 
and thereby reduces the size of m res on the coarse MILC lattices to an acceptable level. 

In Fig. [1] we show our calculated results for m res (as determined from Eq. (HI)) for a 
number of different valence and sea quark values on the coarse MILC lattices. Because the 



2 An analysis of the domain-wall transfer matrix reveals that there is also a contribution to the residual 
quark mass that depends upon L s as p(0)/L s , where p(0) is the density of near-zero eigenmodes [43] . The 
size of both the exponential and power-law contributions to m res have been measured by the RBC and 
UKQCD Collaborations for various domain-wall sea quark actions at multiple lattice spacings [8[. The 
relative importance of the power-law contribution increases with L s and decreases rapidly with the lattice 
spacing. In particular, RBC/UKQCD found this term to be comparable to the exponential contribution 
for L s — 16 at a -1 ~ 1.6 GeV, and 6 times smaller for L s = 16 and a -1 ~ 2.3 GeV. These lattice spacings 
are similar to those used in this work. 
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FIG. 1: Chiral extrapolation of m res on the coarse lattices. The curve shows the extrapola- 
tion/interpolation for points where the domain- wall pion mass is tuned to equal the lightest (taste 
pseudoscalar) staggered pion mass. For comparison, we show the determination of m rea by the 
LHP Collaboration, which uses this tuning [f|. 

mixed-action theory does not have a full QCD (unitary) point, the curve shows an extrapo- 
lation/interpolation to points where the domain- wall pion mass is tuned to equal the lightest 
(taste pseudoscalar) staggered pion mass. For comparison, we show the determination of 
m ies by the LHP Collaboration, which uses this tuning j^]; as one can see, the agreement is 
good. The fit function used was a simple polynomial in the valence and sea quark masses, 
with terms up to fourth order in the valence mass, and linear in the sea mass. The correlated 
fit has a C.L. of around 10%, which is not unreasonable, given that chiral logarithms are 
expected to enter at the same order as the linear terms. We quote a value for rim res in the 
chiral limit on the coarse ensembles of 0.0044(1) (4), where the first error is statistical, and 
the second is a systematic error that comes from varying the functional form and mass range 
used in the fit. In physical units, m res is around 2.7 MeV, close to the size of the physical 
light quark mass. 

In Fig. [2] we compare the values of m res computed on the fine MILC lattices with the 
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FIG. 2: Comparison of rim res on the coarse (a ~ 0.12 fm) and fine (a ~ 0.09 fm) MILC lattices. 
The curves are the interpolation/extrapolation to the LHPC tuning. The value of rim res on the 
fine lattices is approximately three times smaller than it is on the coarse lattices. 



results on the coarse lattices. We observe a significant dependence on the lattice spacing, 
and m res decreases by a factor of three from the coarse to fine lattices. The RBC and 
UKQCD Collaborations see a similar reduction in m res when comparing its value on lattices 
with nearly the same scales as the coarse and fine MILC lattices, but using three flavors 
of domain-wall sea quarks instead of improved staggered sea quarks. We quote a value for 
rim res in the chiral limit on the fine MILC lattices of 0.0016(2), or approximately 1 MeV in 
physical units, where the error is statistical only. 



IV. THE MIXED MESON MASS-SPLITTING 



We have shown that the dominant effect of the residual chiral symmetry breaking is small 
and under control given the size of the residual quark mass, especially on the finer lattices. 
In this section and the following section we demonstrate that discretization effects arising 
due to the interaction of the valence and sea sectors can be understood and controlled using 
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mixed-action chiral perturbation theory (MAXPT). 

MAXPT for Ginsparg- Wilson type quarks on a staggered sea was first considered in 
Ref. [26(, where the one- loop formulas for f w and m 2 were derived. Many additional quan- 



tities have since been calculated to one loop in MAXPT [121 . 1271 . l28l . |29| . It has been 
demonstrated that there is only one low-energy constant particular to the mixed-action chi- 
ral effective theory that appears at one loop. Additionally, only one of the many parameters 
coming from the purely staggered sector enters mixed-action expressions at one loop. Thus 
only two additional parameters enter the mixed-action chiral formulas as compared to those 
for purely domain-wall on domain-wall simulations. These two new parameters are both 
meson mass-splittings, and one has already been determined by the MILC Collaboration 
from spectrum calculations so that it could be used as an input to their chiral fits to more 
complicated quantities [2]. We have calculated the parameter unique to the mixed-action 
case so that we can use this parameter as an input to our own chiral formulas in the same 
manner as the MILC Collaboration. 

In a mixed-action theory one can have mixed mesons made up of one valence and one 
sea quark, in addition to mesons made up of two valence or two sea quarks. We review 
the tree-level mass relations for these mesons as given in Ref. [2a], since they are useful 
in understanding the leading-order lattice-spacing contributions to mixed-action numerical 
simulations. The valence-valence meson mass is given by 

m 2 m , = /idwOu + "V + 2m res ) , (4) 

where we have included the contribution from the residual quark mass. In the limit that 
L s becomes infinite, m ms — ► and we recover the continuum relation. Because of taste- 
symmetry breaking, the staggered theory has 16 pions instead of one. The tree-level relation 
for staggered mesons is [3] 



m as',t = t J Vta S (™> B + my) + A t , (5) 



where the a 2 A t are the taste-splittings of the 16 taste pions. For staggered quarks there 
exists a residual SO(4) taste symmetry at 0(a 2 ), such that there is some degeneracy between 
the 16 pions, and we treat the taste index t as running over the multiplets P, A, T, V, I with 
degeneracies 1, 4, 6, 4, 1. The splitting a 2 A P vanishes because there is an exact (non-singlet) 
lattice axial symmetry. It is convenient to express our formulas in terms of the bare lattice 
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FIG. 3: Determination of the mixed-action parameter A m i x on the coarse and fine MILC lattices. 
The vertical axis is a carefully chosen linear combination of squared meson masses, the left side 
of Eq. (|7j), such that a linear extrapolation of this quantity in the staggered quark mass gives 
the parameter A m j x as the y-intercept. The small diamonds (squares) are the coarse (fine) data. 
For comparison, the large diamonds (squares), show the values of the staggered taste splittings 
measured by MILC on the coarse (fine) lattices. 

masses which enter the lattice calculations, not the renormalized masses. Since the quark 
masses using different actions are renormalized according to different schemes, we absorb 
this scheme dependence into separate coefficients, //dw and // s tag- Finally, the mass of a 
mixed valence-sea meson is 

m 2 vs = fi dw (m v + m res ) + fi staug m s + a 2 A mix , (6) 

where A m i x is the only new constant that appears in the mixed-action theory to next-to- 
leading order. 

We have calculated A mix by computing the mass of a pion made of one domain-wall and 
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one staggered quark. In order to extract A m ; x from our spectrum data, it is convenient to 
construct the following linear combination of squared meson masses: 

m ls - ^ m lv = /Wstag«i s + a 2 A mix . (7) 

We perform a linear fit to the left-hand side of this equation as a function of the staggered 
mass. The result is shown in Figure EJ in the chiral limit this fit gives a 2 A m ; x . Figure [3] also 
compares our results for a 2 A mix with the four taste-breaking staggered meson splittings as 
determined by MILC on both the coarse and fine lattices; we find that A m j x is around the 
same size as the smallest of the MILC splittings. The value of a 2 A mix decreases by roughly a 
factor of two from coarse to fine lattice spacings, which is consistent with the expected 0(a 2 ) 
scaling behavior. We quote a value for r 2 a 2 A mix on the coarse MILC lattices of 0.207(16), 
which is ~ (280 MeV) 2 in physical units, and is consistent with the value found in Ref. [9]. 
On the fine lattices, we find r 2 a 2 A m j x = 0.095(20), which corresponds to ~ (190 MeV) 2 in 
physical units. 



V. UNITARITY VIOLATION AND THE SCALAR CORRELATOR 



Because mixed-action lattice simulations have different discretization errors in the va- 
lence and sea sectors, they are not unitary at non-zero lattice spacing. Such effects are 
expected to vanish in the continuum limit. One must take them into account, however, 
when extracting physical quantities from numerical mixed-action lattice data calculated at 
fixed lattice spacings. This can be done using XPT that is formulated for the specific mixed 
action being used. In the case of domain-wall valence on staggered sea lattice simulations, 
for many quantities of interest such as pseudoscalar meson masses, decay constants, and 
Bk, the violation of unitarity appears mildly in chiral expressions as chiral logarithms of 



the form m 2 ea ln(m 2 val ) 



12 



261 ]. One would like a strong check, however, that MAXPT 



truly describes all unitarity-violating effects that are observable in the lattice data, even 
when the effects are more pronounced. The behavior of the isovector scalar (ao) correlator 
provides such a check. 

Quantities such as the ao correlator and n — n scattering in the 1 = channel are particu- 
larly sensitive to unitarity-violating effects due to the presence of flavor-neutral intermediate 



states 



11 



471 ] . These give rise to diagrams that are disconnected at the quark level, i.e., 
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hairpin diagrams, in which partial quenching effects are more pronounced. In particular, for 
the ao meson, such quark-disconnected diagrams affect the lattice correlator itself. While 
our results in the next subsection are specific to this particular choice of mixed action for 
valence and sea quarks, the issues that we discuss are generic to all mixed-action simulations. 
Similar analyses of the Oq correlator have been done for the two-flavor domain-wall case in 



Ref. 



48 



and for the staggered valence and sea case in Ref. [49]. 



A. The ao in mixed-action XPT 

The isovector scalar is created with the following local operator at the quark level: 

S(x,t) = 3(x, t)u(x, t) . (8) 

In the chiral effective theory that describes the pseudo-Goldstone boson sector, this is rep- 
resented by the point current 

S x (x)=^[^(x)] ud . (9) 

In numerical lattice simulations, the ao correlator is given by C(t) = 
^ x <0|5(x,t)5 t (0,0)|0). The leading contribution to C(t) comes from the propaga- 
tion of an ao meson from time to time t. The propagation of the ao cannot be handled 
within XPT, since in that framework only the light pseudoscalar bosons are dynamical 
degrees of freedom. In XPT the scalar propagator leads to a contact term, but for the 
purposes of our fits, it can be parameterized by an exponential of the form A exp[— m ao t]. 

It was first noticed in the quenched case, however, that the isovector correlator also 

□ 

receives sizeable contributions from two-particle intermediate states [47[ . These two-particle 
"bubble" contributions, unlike the direct exponential term, can be calculated in XPT using 
the expression for the scalar current, Eq. (Q; an example diagram is shown in Figure HI 
Thus the leading expression for the lattice a correlator has two terms: 

C{t) = Ae- m ^ 1 + B{t) + ■■■ (10) 

where the • • • represent excited state contributions. The bubble term B(t) has been calcu- 



lated by Prelovsek using one-loop mixed-action XPT in Ref. The result for 2+1 flavors 



of sea quarks, taking the time direction to be infinite in length but accounting for the finite 
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(til 



do 



du 




FIG. 4: Leading contributions to the scalar current. The first diagram corresponds to the propa- 
gation of an ao meson and the second is one of three "bubble" diagrams. While this figure shows 
a tt and an 77 propagating, there are contributions from tt — tt and K — K intermediate states as 
well. 
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where uof e 

The first important feature of this expression is that there are no free parameters in B(t). 
Both the shape and the normalization are completely predicted by XPT as long as one uses 
point sources for the scalar interpolating field. The meson masses and the coefficient \i can 
be determined from spectrum calculations. The values of mixed-meson splittings A m i x and 
A 1 that appear in m vu , m vs , and mu I , , msi are already known for our choice of valence and 
sea quark actions. As a practical matter, we do not even need the computed residual mass 
as long as we fit B(t) in terms of the meson masses. 

We note that the parameter fi, which relates the quark masses to the squared meson 
masses at leading order in XPT, can be determined in multiple ways. One way is to perform 
a linear fit to m 2 vs. m q . The resulting value of \x = fitme will then be common for all 
valence quark masses at a given lattice spacing. An alternative method is to define \x by the 



3 The form of B{t) shown here is found by integrating the expression in Ref. [n| over the time direct 
well as performing the necessary Fourier transform. 



ion as 
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following ratio: 



2 

m 



/^ratio — 777 ; r, U^J 

2(m„ + m res ) 

which comes from the tree-level XPT expression for the meson masses. This gives rise to a 
different value of /1 for each sea quark ensemble and valence quark mass. We have calculated 
fi using both of these methods, and find that the results are consistent within errors. We 
therefore use /i ra tio for the remainder of our analysis. 

We have considered yet another possibility, which is to take \x in the chiral limit from a fit 
to the pseudoscalar mass data including higher order corrections. By taking ^t rat i (or /i tree ), 
rather than /1 in the chiral limit, one incorporates a subset of higher order corrections to 
the one-loop formulas. This leads to better agreement between XPT and numerical data for 
many purely staggered quantities studied by MILC, including the scalar correlator 0, 
We also find this to be the case in the present work. 

The second important feature of the expression for B(t) is that, unlike in continuum full 
QCD, it receives unphysical contributions from two-pion intermediate states. If the simu- 
lated valence quark masses are sufficiently small, these can dominate the scalar correlator, 
C(t), at large times and make the a mass difficult to determine. Because one of the pions 
in the two-pion bubble is a neutral pion, it has a double pole in its momentum-space propa- 
gator in the isospin limit. The linear- in-t growth factor in the third term of Eq. ( TTTT) comes 
from the resulting double pole, while the minus signs in front of the second and third terms 
come from the hairpin residues. 

Two-pion intermediate states contribute to B(t) whenever the lattice theory is not uni- 
tary; thus they are generically present in any partially-quenched theory. It is only when 
both the masses are equal and the actions are identical in the valence and sea sectors that 
terms 2-4 in Eq. ( TTTT) vanish and B(t) is strictly positive definite. Because full QCD is only 
recovered in the mixed- action theory after taking the continuum limit, two-pion contribu- 
tions to the scalar correlator cannot be removed completely for any choice of mixed-action 
lattice simulation parameters. To illustrate this, we consider the two most "natural" tuning 
choices for the mixed domain-wall valence, staggered sea theory. The first possibility is to 
fix the valence pion mass to be equal to the taste-pseudoscalar sea pion mass. This tuning is 
appealing because the mass of the taste-pseudoscalar pion vanishes in the chiral limit, even 
at finite lattice spacing. This choice, however, increases the unitarity-violating contribu- 
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tions to B{t) because the taste-singlet pion (the only taste in the sea sector which appears) 
is much heavier than the taste-pseudoscalar pion on the coarse MILC lattices. The second 
possibility is to fix the valence pion mass to be equal to the taste-singlet sea pion mass. This 
tuning completely removes the third term in B(t), which has the enhanced (polynomial in 
t) unitarity violation, but still does not eliminate all two-pion contributions. Furthermore, 
this tuning may not be advisable in practice since it would lead to a rather heavy valence 
pion on the coarse MILC lattices. 

Thus, although we could in principle tune the mixed-action theory to remove the enhanced 
unitarity violation in the ao correlator (and in other quantities such as the I = it — it 



scattering phase shift [50J]), the lack of unitarity will always appear to some degree. The 
only way to completely remove behavior such as negative values of the scalar correlator is 
to take the continuum limit, since both rooted-staggered quarks and domain- wall quarks 
should reproduce the same continuum theory as a — > 0. For our purposes, however, a 
large negative bubble contribution to the scalar correlator is actually helpful. As we will 
show in the next section, it allows us to easily see the bubble contribution to our data, and 
therefore to evaluate how well the behavior of the bubble is described by mixed-action chiral 
perturbation theory. 

B. Lattice results for the scalar correlator 

We have calculated the scalar correlator on a subset of the ensembles listed in Table [H 
These are shown in Table [TU along with the valence quark masses and number of configu- 
rations used. Our preliminary calculation of the scalar correlator presented at Lattice 2007 



37| | used point sources and sinks because this is the only case for which MAXPT gives a 
prediction for the bubble term with no unknown parameters. We found, however, that the 
statistical errors using the point source were too large to make a compelling case for the 
agreement between MAXPT and our numerical data, especially on the fine lattices. This 
motivated us to switch to a random-wall source, which simulates many point sources and 
allows us to significantly improve our statistics. The random-wall source has the same nor- 
malization as the point source, up to a factor of the spatial volume, so Eq. ( fill) still applies. 
Fig. [5] shows a comparison of the point source data with the random-wall source data on 
the 0.007/0.05 coarse ensemble with a valence domain-wall mass of 0.01 in lattice units. Al- 
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TABLE II: Data used in the scalar correlator analysis. 



a(fm) ami/ am s am v # configs 

0.09 0.0062/0.0031 0.0062 80 

0.12 0.007/0.05 0.01 184 

0.12 0.007/0.05 0.02 175 

0.12 0.007/0.05 0.03 189 



though both simulations used approximately 200 configurations (206 for the point and 184 
for the random- wall) , the statistical errors in the random- wall data are a dramatic factor of 
5-6 times smaller than those of the point data. 
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FIG. 5: Comparison of the scalar correlator generated with a point source on 206 configurations 
and with a random-wall source on 184 configurations. The data shown is for am„=0.01 on the 
0.007/0.05 coarse ensemble. 

In order to compare our numerical data to MAXPT, we first overlay the prediction for 
the bubble term [Eq. fTTTj) ] on top of the scalar correlator data as a function of time for a 
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FIG. 6: Scalar correlator on the coarse 0.007/0.05 ensemble with three different valence masses. 
Overlaid on the data are the predicted bubble contributions, which should dominate for large time. 



number of different masses and two lattice spacings. Fig. [6] shows three different valence 
masses on the 0.007/0.05 coarse ensemble, while Fig. [7] shows two similar valence masses 
but on different lattice spacings. The prediction for the bubble term in Eq. (iTTj) requires 
as inputs three parameters of the chiral effective theory: A mix , which we have calculated in 
Section IV; A/, which we take from Ref. 3|; and /i, which we obtain from our mixed-action 
pseudoscalar data using Eq. (fl2l) . The bubble term is expected to dominate the scalar 
correlator at large times, and in this region the agreement between the data and the XPT 
prediction is good. For sufficiently small times, the direct term Ae~ ma Q t dominates, and 
we would not expect the bubble to match the data in this region. Note that the trends in 
the data are qualitatively what one would expect based on Eq. (II ip ; Fig. O shows that, as 
the valence quark mass is made lighter, the negative contribution to the bubble increases 
in magnitude at fixed sea quark mass, and the correlator C(t) becomes negative for much 
of the time range for the smallest mass. Fig. [7] illustrates the lattice spacing dependence 
of the scalar correlator and the predicted bubble term. The negative contribution to the 
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FIG. 7: Scalar correlator on the coarse 0.007/0.05 ensemble for am„=0.01 and on the fine 
0.0062/0.031 ensemble for om„=0.0062. As in Fig. [6l we include the prediction for the bubble 
expression for comparison. Note that we have rescaled the time axis for the results on the fine 
lattice so that we can compare the two different lattice spacings. The y-axis is dimensionless. 



bubble is much smaller on the fine lattice than on the coarse even though the valence- valence 
and taste pseudoscalar sea-sea pion masses are roughly comparable in physical units; this is 
because the mass-splittings a 2 A mix and a 2 A/ are smaller. 

If we ignore lattice artifacts and use the continuum form of the bubble expression (by 
setting a 2 A/ = a 2 A mix = 0), we find that the bubble contribution never becomes negative for 
the masses used in this work, which is in clear disagreement with the data. The continuum 
partially-quenched bubble is plotted in Fig.[H]for the parameters corresponding to our lightest 
valence mass on the coarse lattice, along with the actual numerical data and the bubble 
prediction including the unitarity-violating lattice artifacts. The continuum curve does not 
agree with the data, while the curve that includes the mass-splittings is in good agreement 
for large times. This further emphasizes the need to include lattice artifacts such as the 
mass-splittings in the chiral effective theory in order to describe correctly the lattice data. 
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FIG. 8: Comparison of the bubble prediction for the ami/am s = 0.007/0.05, am v = 0.01 mixed- 
action data with and without the mass-splittings . 

As a consistency check of the chiral parameter /i that is used in the bubble prediction, we 
fit the various data sets with \i free at large times where we can ignore the direct exponential 
contribution to the scalar correlator. The results of these fits are shown in Table IIHI We 
find that the values for /i determined from the fits are consistent with the values of /i ra tio 
used in the prediction of the bubble term at the ~ 1.5<r level. For the valence mass 0.03 on 
the coarse ensemble, the bubble term is negligible compared to the statistical errors on the 
data, and this test is not possible. Fig. [9] shows a plot of the scalar correlator for the lightest 
valence mass on the coarse ensemble with a fit to the bubble expression at large times. 

Having confirmed that the results for [i from fits to the scalar correlator are consistent 
with the values of // ra tio determined from pseudoscalar data, we next perform fits to the full 
expression for C (t) in Eq. ffTUl) . Table [IV] shows the results of the different fits for each of 
the ensembles, including the bubble term with \i set to /^atio plus a ground state exponential 
term. In order to describe the fine lattice data at small times, we find that we must include 
an additional exponential term to model excited state contributions. We can get reasonable 
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ami am v a/i free a/i rat io (*min, ^max) X 2 /d.o.f. 



0.0062 0.0062 0.86(14) 1.070 (12,25) 1.24 
0.007 0.01 1.81(11) 1.629 (6,20) 0.98 
0.007 0.02 1.42(33) 1.541 (6,20) 1.46 

TABLE III: Fits to B(t) at large times leaving fi as a free parameter. For comparison we also show 
the value of /i ra tio a s determined by Eq. (fl~2l) . 
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FIG. 9: Fit of the ami/am s = 0.007/0.05, am v = 0.01 data at large times to the bubble expression, 
B(t), keeping ^ as a free parameter. The fit range was 6 to 20. 



fits to all masses shown, although we point out that for the heavier mass of m v = 0.03 on 
the coarse lattice, the bubble contribution has a very small effect [for t > 4, B(t) ~ 0], 
and similar \ 2 results can be obtained while omitting it. For the lighter masses, including 
the bubble is absolutely necessary. Figs. [TO] and [H] show fits of the coarse data to the 
full expression for C(t) including only the ground-state exponential plus the bubble term. 
Fig. [TD] shows the results for all three valence masses on the coarse ensemble, while Fig. [TT1 
shows only the lightest valence mass, but Fig. [11] includes the error band for the resulting 
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fit. Given our current statistics, we can successfully fit all of the coarse data starting at a 
minimum time of three or four while including only a single exponential term. This is not 
the case, however, for our fine lattice data which has both a smaller bubble contribution and 
a finer resolution in the time direction. When we include an excited state term in the fit 
to the fine lattice data, we are able to obtain an acceptable x 2 /d.o.f. over a much larger fit 
range, and we reduce the statistical errors in the fitted parameters, as illustrated by Fig. [T2l 
The different fits of the fine lattice data yield consistent values for the ground state meson 
mass, and we find am ao « 0.38(5)(3), where the first error is statistical and the second is 
a systematic error which reflects the differing central values between the single and double 
exponential fits. Given that a~ x ~ 2.3 GeV on the fine lattice, we note that this result is 
close to the experimentally measured a® meson mass of ~ 980 MeV. We emphasize, however, 
that this result is for unphysical quark masses at a nonzero lattice spacing, and that the 
systematic error does not include other important sources of uncertainty such as quark mass 
dependence and discretization effects. 
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FIG. 10: Fits to all three valence quark masses on the 0.007/0.05 coarse ensemble. The fit function 
is given in Eq. (|10p . and contains a single exponential plus the bubble contribution with /i — /^ratio- 
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FIG. 11: Fit to the ami/am s =0. 007/0. 05, 077^=0. 01 data with error band displayed. Note that 
this is the same fit as the solid black curve in Fig. [TUJ 



#Exp 


ami o- m v 


(^minj tmaxj 


X 2 /d.o.f. 
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0.0062 0.0062 


(7,24) 


1.10 
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0.0062 0.0062 


(2,21) 


0.87 


1 
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1.04 
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(4,20) 


0.89 
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0.57 



TABLE IV: Fits to the mixed-action scalar data using the expression for C{t) given in Eq. (]10p . 
keeping \i fixed to the ratio value defined in Eq. (|12f) . For the fine lattice data only, we show results 
for both one and two exponential terms in the fit function. We do not show the resulting masses 
and amplitudes in this table because our purpose is not to calculate the ao mass, but rather to 
demonstrate that we can successfully fit the data with a physically motivated fitting function. 
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FIG. 12: Comparison of fits to the fine scalar data with a single (grey checked) or double (cyan solid) 
exponential term. The single exponential fit required a t m i n of 7 in order to obtain a reasonable 
X 2 /d.o.f.; the double exponential fit allowed a t m i n of 2. 



Although it has been suggested in Ref. [Ill that the scalar correlator could be used to 
determine A mix , we have found this to be very difficult with the parameters and level of 
statistics used in the present work. We note that varying A m j x within the statistical errors 
from Sec. IIVI or even by a factor of two, does not have a significant effect on the numerical 
value of the bubble term for our values of the simulation parameters. Consequently, since 
the bubble term is not very sensitive to this variation, we would not expect A m i x to be 
well-constrained by our scalar data, and fits leaving A m i x as a free parameter confirm this 
expectation. Finally, we observe that the bubble term is much more sensitive to Aj, and the 
consistency of our results with MAXPT provides a good check of the MILC determination 
of A/. 

We conclude this section by observing that our data is quite well described by the mixed- 
action XPT formula calculated by Prelovsek. This highly non-trivial test indicates that 
mixed-action XPT can describe the largest source of unitarity violations that appear in our 
simulations. Additionally, our data shows that the degree of unitarity violation decreases 
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with the lattice spacing, as expected. We note, however, that other sources of unitarity 
violations are expected to be present at some level, although we cannot currently resolve 
them in our data. These are due to short- distance effects and therefore cannot be described 
by XPT. The first of these effects is due to the fact that domain-wall fermions do not have 
a positive-definite transfer matrix 411 ] . This leads to the appearance of opposite parity 
states in two-point correlators, which are similar to typical excited states, except that they 
oscillate in sign from one time slice to the next. Evidence of these oscillating states using 
the same mixed-action scheme as this work was presented in Ref . 5l| , where the effect was 
shown to be small in the pion two-point correlator. We note that oscillating states also 
contribute in the staggered valence on staggered sea case, and the analysis of the scalar 



correlator in Ref. 



49j required such terms in order to describe the numerical data. We do 



not see evidence of oscillating states, however, in our own scalar correlator data. The second 
of the unitarity-violating effects is the presence of enhanced zero-mode contributions in the 
(partially) quenched theory; this was studied in Ref. [451 ] . These effects are known to be 
suppressed by the fermion determinant, large spatial volumes, and large time separations. 
Although such effects must be present in partially-quenched domain-wall correlators at some 
level, we see no evidence for them in this work. As we have shown, our numerical mixed- 
action scalar data is quantitatively consistent with the MAXPT bubble prediction plus a 
small number of non-oscillating exponential terms. 



VI. CONCLUSIONS 



Mixed-action lattice simulations using domain-wall valence quarks and improved stag- 
gered sea quarks are well-suited to the calculation of weak matrix elements such as B^- 
Discretization effects in the mixed-action theory arise from both the valence and sea sectors, 
as well as the interaction between the two. Although such effects vanish in the continuum 
limit, in practice they must be removed from quantities calculated at finite lattice spac- 
ing with the aid of extrapolation formulae calculated in mixed-action Chiral Perturbation 
Theory. 

In this work we have studied the size of discretization effects in the LHP Collaboration 
mixed-action scheme which uses the MILC Asqtad-improved staggered lattice configura- 
tions. We have calculated the values of the residual quark mass, m res , and the mixed 
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valence-sea meson mass-splitting, A m ; x , and have found them to be comparable in size to 
or smaller than analogous parameters in numerical simulations by the RBC, UKQCD, and 
MILC Collaborations. Thus we conclude that the size of generic discretization effects in the 
mixed-action theory are small enough to allow precise determinations of continuum weak 
matrix elements. We have also performed a strong check of the ability of MAXPT to ac- 
curately describe unitarity-violating discretization effects in the isovector scalar correlator. 
We find that the MAXPT prediction for the two-particle intermediate state (bubble) contri- 
bution to the scalar correlator is in good quantitative agreement with the numerical lattice 
data. Thus we conclude that MAXPT correctly describes the dominant unitarity-violating 
contributions to mixed-action lattice simulations. 

In the case of the scalar meson, the unitarity-violating discretization effects are much 
larger than the meson ground-state signal for many of our quark mass values, so any un- 
certainty in the bubble contribution translates directly into an uncertainty in the extracted 
scalar mass. One may therefore argue that mixed-action (or any partially-quenched) sim- 
ulations cannot be used for a clean determination of the scalar meson mass. It was shown 
in Ref. 4B 



however, that for the case of Nf = 2 domain-wall fermions, use of partially- 
quenched data allows a more precise extraction of the mass than with full QCD points 
alone, despite contamination from unphysical intermediate states. In fact, it is interesting 
to note that, because MAXPT predicts the observed unitarity violations so well, one can use 
the generic framework of partially-quenched lattice XPT to aid in choosing the best actions 
and parameters for determining the a® mass. We have not done this because our primary 
goal is the calculation of weak matrix elements. Fortunately, in the case of Bk, MAXPT 
predicts that non-analytic unitarity-violating errors should contribute less than a percent 
on the coarse MILC lattices. This fact, in conjunction with our successful analysis of the 
scalar correlator, substantiates the claim that unitarity-violating effects in mixed-action lat- 
tice simulations can be accounted for and removed to recover precise continuum values for 
weak matrix elements. 
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